Doping Dependence of Bilayer Resonant Spin Excitations in (Y, Ca)Ba 2 Cu 3 06+ x 
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Resonant magnetic modes with odd and even symmetries were studied by inelastic neutron scat- 
tering experiments in the bilayer high-T c superconductor Yi_ x Ca x Ba2Cu3C>6+y over a wide doping 
range. The threshold of the spin excitation continuum in the superconducting state, deduced from 
the energies and spectral weights of both modes, is compared with the superconducting d-wave 
gap, measured on the same samples by electronic Raman scattering in the B\ g symmetry. Above 
a critical doping level of 5 ~ 0.19, both mode energies and the continuum threshold coincide. We 
find a simple scaling relationship between the characteristic energies and spectral weights of both 
modes, which indicates that the resonant modes are bound states in the superconducting energy 
gap, as predicted by the spin-exciton model of the resonant mode. 



In high-T c copper oxides superconductors, inelastic 
neutron scattering (INS) experiments have shown that 
the superconducting (SC) phase exhibits an unusual 
spin triplet excitation, the so-called magnetic resonance 
mode yj. This excitation has been observed in sev- 
eral families of copper oxides with SC critical temper- 
atures T c -90 K: Tl 2 Ba 2 Cu0 6+x , Bi 2 Sr 2 CaCu 2 8+x , 
and YBa 2 Cu306+ x (YBC06+:c)- Two kinds of theoret- 
ical models attempt to describe this mode. Within 
an itinerant-electron picture, it is described as a spin- 
exciton, that is, it corresponds to a spin-triplet bound 
states stabilized by the electron-electron interaction in 
the plane below the electron- hole (e-h) spin- flip contin- 
uum, which is gapped in the SC state |2J,|3|]. As the mode 
is an intrinsic feature of the SC state, it disappears above 
T c . In other approaches 0,0, |(J, the collective spin ex- 
citations pre-exist in the normal state, where they are 
damped by scattering from charge excitations. In the SC 
state, the damping is suppressed, and the mode sharp- 
ens if its energy is below the threshold of the gapped e-h 
continuum. The charge carriers interacting with the spin 
excitations can be either uniformly distributed Q| or seg- 
regated into quasi-one-dimensional stripes that separate 
nearly insulating antiferromagnetic domains 0, ||J . Al- 
though the quantum numbers characterizing the mode 
are identical in both of these competing approaches, the 
underlying physical pictures are quite different. It is 
hence important to establish criteria that allow quan- 
titative experimental tests of these models. 

The bilayer structure of YBCO offers an interesting 
opportunity in this regard. While most calculations have 
been carried out for a single Cu0 2 plane per unit cell 
0, 0, both theoretical approaches predict two distinct 
magnetic resonant modes in a bilayer system. They are 
characterized by even (e) or odd (o) symmetries, re- 
spectively, with respect to the exchange interaction be- 
tween the two layers in a bilayer unit. Recent INS mea- 
surements have shown that the bilayer system YBCO 
indeed exhibits two magnetic resonant modes, both of 



which are apparent as a strong enhancement of the mag- 
netic intensity at the planar antiferromagnetic wave vec- 
tor qAF = (tt, tt) in the SC state Their intensi- 
ties display an order-parameter-like temperature depen- 
dence and vanish at T c , without any significant shift 
of their characteristic energies. A crucial difference be- 
tween both theoretical approaches concerns the energy- 
integrated spectral weight (SW) of the magnetic modes. 
In spin exciton models, this quantity is a function of the 
binding energy of the modes j2j, j3|(that is, the difference 
of the mode energy and the e-h continuum threshold), 
whereas in models based on pre-existing spin excitations 
it is not obviously related the location of the mode with 
respect to the e-h continuum. 

Here we report an INS study of the odd and even res- 
onant spin excitations at qAF spanning a wide doping 
range from the underdoped to the overdoped regimes in 
YBCO. The threshold for e-h excitations, 2A max , was 
measured by electronic Raman scattering (ERS) in the 
Big channel 9] on the same samples used in the neutron 
experiments. This eliminates systematic errors invariably 
associated with comparisons of measurements on samples 
from different origins. We found that the resonant spin 
excitations are always located below 2A max . Assuming 
that the spin excitations are bound states in the energy 
gap of the e-h continuum, its threshold at qAF, k> c , is 
estimated using the measured energies, E°' e , and spec- 
tral weights, W°' c , of the odd and even resonance peaks. 
The threshold energy lu c estimated in this way is found 
to be below the independently determined gap 2A max , 
so that the analysis is self-consistent. We also observed 
a systematic scaling relation between W°'° and the re- 
duced binding energy (w c -E° ,c )/w c . As such a relation is 
only expected for collective bound states of the e-h con- 
tinuum, our results support the spin-exciton description 
of the resonant mode. Further, our study reveals that 
2A max and the total SW of the resonant spin excitations 
at qAF drop precipitously in the overdoped regime close 
to the hole doping level 6 ~0.19. Previous work identified 
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of this do ping lev el as the end point of the pseudo-gap 
phase mHlHli- 

In addition to previous INS experiments on a 
slightly overdoped sample (OD85, Y .gCao.iBa2Cu307, 
Ref. 0) and a slightly underdoped sample (UD89), 
YBa 2 Cu30 6 . S5 , Ref. 0), INS studies of bilayer excita- 
tions were performed on two YBCO samples. The first 
sample is overdoped (OD75), Yo.ssCao.isE^CuaOy with 
T c = 75 K0 Il4j . Following our previous work with 
sample OD85, about 50 single crystals, obtained by a 
top-seeded solution growth method 0], were co-aligned 
on aluminum plates and fixed with glue or Al screws. The 
second sample (UD63), YBa 2 Cu 3 6 . 6 with T° nset = 63 
K, is underdoped and consists of 180 detwinned squared- 
shape single crystals 0|. The mosaic spread of these 
crystals arrays exceeds 1.4°. INS measurements have 
been performed on the thermal triple axis spectrome- 
ter IN8 at the Institute Laue Langevin (Grenoble). The 
experimental setup and scattering plane are the same as 
Ref. except that we utilized a PG(002) monochroma- 
tor. A fixed final neutron energy of 35 meV has been 
used yielding a typical energy resolution of 7-8 meV. 

The spin susceptibility of a bilayer system reads: 

x(Q,w)= sin 2 (7T2:L)x (q,w) +cos 2 (7rzi)x e (q,w) (1) 

where Q = (H, K, L) is the full wave vector and q = 
(H,K) is the planar wave vector in CuC>2 plane. z=0.28 
stands for the reduced distance between the planes of the 
bilayer. The L-component of the wave vector can thus 
be used to select either the odd or the even channel ||: 
the magnetic structure factor is maximum for 5.2 in 
the odd channel and for 3.4 or 7 in the even. Follow- 
ing the standard terminology H, 0, El , we define the 
magnetic resonant mode as the enhancement of the imag- 
inary part of the magnetic susceptibility, A7mx(qAF, w) 
in the SC state, which is derived from the difference be- 
tween measurements in the SC state at 10 K and the 
normal state above T c . The imaginary part of the dy- 
namical magnetic susceptibility is also calibrated in abso- 
lute units (/i|.eV _1 /f.u.) following a standard procedure 
with a reference phonon at 42.5 meV [rR lisj. 

Fig. provides evidence of the existence of odd and 
even resonant magnetic peaks in the INS spectra of the 
highly overdoped sample OD75. These results differ sig- 
nificantly from previous reports at lower doping levels 
0] , because the energies of both mode are very sim- 
ilar: E°=34meV and E^=35meV. Further, both modes 
also display similar amplitudes. In order to determine 
the amplitudes I°' e of both modes, the L-dependence of 
the magnetic intensity at 34 meV (Fig. \l\b) was fitted 
with Eq. n This yields a ratio I c /I° ~ 0.4, compared 
to 1/3 for sample OD85 and 1/6 for YBC0 6 . 85 @ 
However, as for lower doping levels 0@], the tempera- 
ture dependence of the intensities of both modes exhibits 
a marked change at the SC temperature (Fig. ^c), in- 




-400 tf 

20 30 40 50 

Energy (meV) 




34 meV J 3050 d ) 

— I — R A - 



too/ 

(-0.5,-0.5,L) 

2500 

UD63 10K-70K 

, ■= a- 2000 
— n — L=5.4 



20 40 60 80 100 

Temperature (K) 




30 40 50 60 70 80 

Energy (meV) 



FIG. 1: (color online) a) Differences between INS energy 
scans measured at the antiferromagnetic wave vector be- 
low and below T c in Yo.85Cao.i5Ba2Cu307 (OD75). Open 
and closed symbols correspond to the odd and even chan- 
nels, respectively. The data have been corrected by the 
anisotropic Cu magnetic form factor. The stars indicate the 
measured reference level of the magnetic scattering deter- 
mined by constant-energy scans, b) L-dependence at E=34 
meV for OD75. The lines describe the intensity of each mode 
following Eq. multiplied by the Cu magnetic form factor 
(upper line), c) Temperature dependence of the neutron in- 
tensity measured at characteristic energies of the even and 
odd modes for OD75. d) Same as a) but for sample YBC06.6 
(UD63). 



dicating that the phenomena have a similar origin. In 
the strongly underdoped sample UD63 (Fig. \J\d), the 
two mode energies and their intensities are very different: 
E°=37 meV, E£=55 meV, and I c /I° ~ 1/20. As already 
observed in the OD85 sample |2|, the odd mode exhibits 
a broader energy line shape in OD75, AE~ 12 meV, than 
the resolution-limited energy profile observed in under- 
doped samples [16]. The energy width of the even mode, 
AE^ 15 meV, is always broader than the resolution, for 
all doping level. 

In order to assess whether or not these resonant modes 
should be regarded as bound states in the electronic con- 
tinuum, one has to determine the amplitude of the SC 
gap. The good surface quality of the single crystals re- 
moved from the crystals array OD85 and OD75 enables 
ERS experiments. The ERS technique probes electronic 
excitations in selected areas of the Brillouin zone. In 
the Big channel, the antinodal areas around (±7r,0) and 
(0,±7r) are probed, and a peak appears in the electronic 
spectrum in the SC state at twice the energy of the 
maximum e?-wave SC gap, 2A max j^]. These measure- 
ments were carried out with a triple-grating spectrome- 
ter in quasi-back-scattering geometry. The crystals were 
mounted on the cold finger of an He circulation cryostat. 
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FIG. 2: (color online) Electronic Raman spectra in the 
Bigchannel. The difference between spectra measured at 15 K 
in the SC state and at 90 K in the normal state is displayed. In 
addition to sharp features due to phonon renormalization, the 
difference indicates a broad peak due to the enhancement of 
the electronic response in the SC state. The electronic peaks 
were fitted with Gaussian profiles. The inset shows the peak 
energy of the 5 samples we have studied (see Fig. caption 
for the determination of the doping levels). 

The 514 nm excitation line of a Ar + /Kr + mixed-gas ion 
laser was used. 

Figure shows the Bi g peaks in the ERS spectra of 
single crystals extracted from samples OD85 and OD75. 
The Big peak is observed at 69.3 meV in the first sam- 
ple and shifts to 35.5 meV for the second, whereas T c 
is only reduced from 85.5 K to 75 K. The same results 
were obtained on other crystals selected randomly from 
the OD85 and OD75 arrays[2l]. Further experiments on 
samples with different T c also confirm that the charac- 
teristic energy of the Bi g peak drops steeply by almost 
a factor two with a minor reduction of T r . This phe- 
nomenon was previously reported in Refs. 19ll20|, and it 
is now confirmed directly on the same samples used for 
the INS experiments. Systematic errors associated with 
experiments on samples of different origins could thus 
be avoided. As shown in the inset of Fig. [21 the drastic 
decrease of the Bi g peak occurs close to a critical hole 
doping level S c ~ 0.19. 

In Fig. 0a, we report the hole doping dependence (i) of 
the of 2A m ax determined by ERS from the present study 
and Refs. 0, |2(J, and (ii) the odd and even resonant 
energies obtained from INS data. For all samples, the 
even resonance appears at higher energy and with less 
intensity than the odd resonance. While E° evolves as 
a function of hole doping as ~ 5 keT c (Fig. 0a), E° 
does not scale with T c : it remains almost constant in the 
underdoped regime and begins to decrease upon entering 
the overdoped regime. With increasing hole doping, the 
distance in energy between the magnetic modes decreases 
from 17 meV in UD63 to 1 meV in OD75. Owing to the 
differences between the energy lineshapes of the odd and 
even modes, the most meaningful comparison is based 
on their energy-integrated spectral weights: W° ,e = J Q 
dwA7mx 0)e (q, u>). The ratio and the sum of the energy- 
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FIG. 3: (color online) Doping dependence of: a) E°, E°, uj c 
and 2A max (ERS). Other ERS data were obtained from [2Q] 
(A), in 19] (o), in the present study (v)- b) Ratio and sum 
of the energy-integrated SW. The hole doping is estimated 
from thephenomenological relation T c (<5) = T° pt (1 - 82.6(5- 
<5o P t) 2 ) H, Sopt is generally estimated to be 0.16 0. T° pt 
is 93K for pure YBCO. Ca substitution reduces T c recording 
to T° pt = 90.6 — 39.2xca obtained for our single crystals. 

integrated SWs of both modes are shown in Fig. 0b. 
Obviously, /W° increases with increasing doping. 

The dramatic change of the SC gap measured by ERS 
(Fig. 0a) separates two distinct regimes on both sides 
of S c : in regime I, encompassing the underdoped sam- 
ples up to OD85, one obtains the following hierarchy: 
5k B T c - E° < E° < 2A max . In regime II (OD75), 
all energies collapse to 5kBT c ~ E° ~ E° ~ 2A max . 
The coincidence of the resonant mode energy and the 
SC energy gap in the overdoped range was suggested 
by prior tunneling experiments |22j . As the information 
about the mode energy obtained from tunneling is quite 
indirect, the interpretation of these data has remained 
controversial. The total SW of the resonant modes, 
W° + Wf:, declines precipitously around a similar doping 
level (Fig. 0b). It is interesting to note that the doping 
level S c ~ 0.19 separating these two regimes coincides 
with the hole concentration where previous work had un- 
covered rapid doping-induced modifications of the phys- 
ical properties such as resistivity, specific heat, Knight 
shift, and muon spin relaxation rate [lfj, [Tl|, [l^] • Accord- 
ing to |23| , this point in the phase diagram is a quantum 
critical point associated with the doping-induced disap- 
pearance of a putative order parameter characterizing the 
pseudo-gap phase. Therefore, it is tempting to relate the 
anomalous hole doping dependencies shown in Fig. a to 
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FIG. 4: (color online) a) Resonant modes SW vs the reduced 
binding energy. Sketch of two approaches for the resonant 
modes: b) bound states of the electronic continuum, c) pre- 
existing modes whose intensities are enhanced by the disap- 
pearance of the fermionic damping below uj c . 



itative differences between the excitation spectra in the 
SC and normal states at energies comparable to the SC 
gap pjj . In overdoped samples, normal-state excitations 
appear to be heavily overdamped and have thus far not 
been clearly identified. It is therefore reasonable to con- 
sider the excitations in the SC state separately, as we 
have done. 

In summary, odd and even symmetry resonant mag- 
netic excitations have been observed in bilayer YBCO 
over a wide doping range. A detailed analysis of the 
modes energies and spectral weights show that the res- 
onant mode in the high-T c cuprates arises from bound 
states in the superconducting energy gap. Our results 
have no natural explanation in models where the resonant 
excitations are associated with excitations pre-existing in 
the normal state. 



the pseudo-gap phenomenon. 

In Fig. [3Ja, one also notices that the resonance en- 
ergies are systematically lower than twice the SC gap. 
This indicates that the resonant modes are below the 
continuum threshold at qAF, w c , which can be estimated 
as 1.8A max using the typical Fermi surface topology ob- 
served in cuprates |2^]- As discussed in Refs. 0,O|, an- 
other estimate of u> c can be made within the spin-exciton 
model where W°' e oc (uj c — E°' e ) /oj c . lu c is simply deduced 
by the ratio of the SWs of the odd and even modes 0, @ • 
For the four samples, uj c extracted in this way from INS 
(Fig. 01a) agrees with the value deduced from ERS data. 

We now discuss these observations in the light of 
the various models attempting to describe the resonant 
modes 0- C3- d 0- E3l As sketched in Fig. one can en- 
visage two situations: either the modes are bound states 
in the energy gap of the electronic continuum (Fig. QJb), 
as assumed in the spin-exciton model |2|. L|- or they al- 
ready exist above T c 0, |f| Q and their SWs are only 
enhanced in the SC state (Fig.^Jc). In the former case, 
the SWs of the modes are proportional to the binding 
energies, while in the latter case there is no specific rela- 
tion between the SW and the distance between the mode 
energy and u> c . The SWs of the modes, W°' e , are pre- 
sented in Fig. 0] function of the reduced binding 
energy, i.e. (uj c — E°< e )/u> c . Clearly, these two quantities 
are proportional to each other with a constant scaling 
factor over a wide doping level. This is the expected be- 
havior for collective bound states in the SC energy gap 
where the scaling factor is mostly controlled by doping 
independent microscopic parameters. In the pre-existing 
mode picture, on the other hand, the SW ratio /W° 
would be also determined by the spin dynamics in the 
normal state. To be sure, well-defined magnetic modes 
are observed in the underdoped samples above T c [l7| . 
but we recently found that there are a number of qual- 
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